Spatial and temporal tuning in void models for acceleration 
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There has been considerable interest in recent years in cosmological models in which we inhabit 
a very large, underdense void as an alternative to dark energy. A longstanding objection to this 
proposal is that observations limit our position to be very close to the void centre. By selecting 
from a family of void profiles that fit supernova luminosity data, we carefully determine how far 
from the centre we could be. To do so, we use the observed dipole component of the cosmic 
microwave background, as well as an additional stochastic peculiar velocity arising from primordial 
perturbations. We find that we are constrained to live within 80 Mpc of the centre of a void — a 
somewhat weaker constraint than found in previous studies, but nevertheless a strong violation of 
the Copernican principle. By considering how such a Gpc-scale void would appear on the microwave 
sky, we also show that there can be a maximum of one of these voids within our Hubble radius. 
Hence, the constraint on our position corresponds to a fraction of the Hubble volume of order 10 -8 . 
Finally, we use the fact that void models only look temporarily similar to a cosmological-constant- 
dominated universe to argue that these models are not free of temporal fine-tuning. 

PACS numbers: 98.80.Es, 95.36,+x, 98.65.Dx 



I. INTRODUCTION 

Modern observational data suggest that the standard 
model of cosmology (SMC; see, e.g., [1]) now rests on firm 
foundations. Two of its main postulates are the Coperni- 
can principle and the existence of dark energy that closely 
resembles a cosmological constant. Nevertheless, inho- 
mogeneous models that propose to do away with both of 
these ideas have been studied as alternatives to the SMC 
for over a decade (see [2-4] and references within for re- 
cent studies, or [5] for a more general review). These 
proposals aim to dispose of the perceived inelegance of 
dark energy by appealing to the fact that we are usually 
restricted to observing the Universe along our past light 
cone. In particular, an increasing expansion rate in time 
[the standard explanation for the redshift-luminosity dis- 
tance relation of recent surveys of Type la supernovae 
(SNe)] could also be interpreted as an expansion rate 
that increases towards us spatially if the assumption of 
large-scale homogeneity were to be abandoned. Indeed, 
supposing that we live in a large (^Gpc scale) under- 
dense region that is isotropic about us, it is possible to 
obtain an exact match to the redshift-luminosity distance 
relation of the SMC (see, e.g., [6]). 

While this is an encouraging sign, there are several ob- 
stacles to a wholesale adoption of such a "void" model 
as a solution to the mysteries surrounding dark energy. 
Voids of sufficient size occur with negligible probability in 
our current picture of structure formation [7], and must 
be nearly spherical to accord with the observed isotropy 
in the cosmic microwave background (CMB). Further 
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constraints abound, in the form of baryon acoustic os- 
cillations [8], local estimates of the Hubble rate [9], the 
kinematic Sunyaev-Zeldovich effect [10], spectral distor- 
tions in the CMB [11], and possible future measurements 
of cosmic shear [12, 13]. The information contained in 
the full CMB power spectrum has been found to place 
very strong constraints on void models [2, 12]. 

A different kind of constraint to consider is the max- 
imum distance an observer could be from the centre of 
a void to match current observations. A rough guideline 
can be obtained by calculating the anisotropic redshift- 
luminosity distance relation that would result from an 
off-centre position [14, 15]. However, it is well-known [16- 
19] that an off-centre observer would observe a prominent 
dipole component of the CMB, since photons traveling 
from the last scattering surface (LSS) would experience 
direction-dependent degrees of redshifting. 

Previous studies [3, 14, 20, 21] have asked how far 
from the centre an observer could live before detecting 
a dipole exceeding the currently measured value. Our 
primary aim is to perform a thorough calculation of this 
distance, distinct in several ways from earlier estimates. 
First, we examine a broad subset of the space of possible 
void profiles, building upon the framework established in 
Refs. [2, 12]. Second, we use the CMB dipole measured 
for the Local Group («lg) m obtaining our constraint, 
rather than, say, the dipole for the solar system, for rea- 
sons discussed below. Finally, we add an extra stochastic 
peculiar velocity component to the dipole induced by an 
off-centre position, arising from the perturbed velocity 
field created by density fluctuations at early times. 

We begin in Sec. II by reviewing the formalism used 
in our investigation. In Sec. Ill A, we present compu- 
tations demonstrating that no other large voids can be 
present within our Hubble volume. Thereafter, we de- 
scribe the computational methods employed, and state 
our results for both radial (Sec. IIIB) and volume (Sec. 
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III C) constraints. We also find in Sec. IV that the 
redshift-luminosity distance relation for supernovae im- 
poses some fine-tuning of the time at which a void is 
observed. We conclude in Sec. V with a summary of our 
work. 



II. FORMALISM 
A. Induced dipole 

Following standard practice, we consider the void to be 
described by a Lemaitre-Tolman-Bondi (LTB) spacetime 
[22-24]: an exact, spherically-symmetric dust solution 
to Einstein's equations. Details of this solution and our 
numerical implementation of it are provided in Ref. [2]. 
We first wish to calculate the CMB dipole induced by 
the void at an off-centre position, and then incorporate 
an additional stochastic peculiar velocity. The induced 
dipole is found using the method described in [2] . Briefly, 
the redshifts for two radial geodesies are considered, one 
incoming from the LSS (at z- m , set to 1000) to the ob- 
server, and the other outgoing through the void centre 
in the opposite direction (z out ). The outgoing ray is con- 
strained to originate at the same proper time (and hence 
the same matter density) as the incoming ray. 

Setting c = 1, the resulting dipole velocity v r can be 
found via 

v r = WT 7 ' W 

Z T" 2>otlt T" Z[ n 

This determines the temperature anisotropy through 

^=v T co S 6, (2) 
or, equivalently, the dipole amplitude via 




In general we expect higher multipoles to appear at 
higher order in w r , but calculations involving the full non- 
radial geodesies confirm that the dipole contribution is 
dominant at the distances from the centre that we con- 
sider [20]. Numerically, our technique is found to agree 
with the analytical approach given in [21] to better than 
one percent. 

Note that our calculation of v r is equivalent to a cal- 
culation of the relative tilt between the comoving matter 
and radiation frames at the off-centre observation point. 
Importantly, the LTB solution assumes a purely dust 
source, and hence cannot properly take into account the 
effect of radiation on the background evolution. However, 
at late times the radiation density is much smaller than 
the matter density, and hence it is a good approximation 
to treat the radiation as a test field on the LTB space- 
time, as we have done. At early times (or sufficiently far 



outside the void), the void spacetime approaches homoge- 
neous Friedmann-Lemaitre-Robertson- Walker (FLRW) , 
since we only consider growing mode void profiles. Thus, 
again, in this regime we do not require the full relativis- 
tic solution including radiation and matter in spherical 
symmetry. 

B. Stochastic peculiar velocity 

Due to the complexity of describing perturbations 
on general LTB backgrounds, directly evolving the per- 
turbed velocity to find the extra stochastic component of 
the dipole is not currently feasible. However, following 
the approach of [2], we note that close to the void cen- 
tre, shear is small enough that perturbations of matter 
density and expansion evolve essentially as in an open 
FLRW model. Thus, for the small scales making the 
dominant contribution to the dipole, we are able to ap- 
ply the tools of Newtonian perturbation theory on FLRW 
backgrounds: for a single spatial component Vi of the ran- 
dom peculiar velocity, the power spectrum reads [25] 

= ^£^M*)> (4) 

where / ~ f2~; ^ s a factor accounting for the suppres- 
sion of growth of velocity fluctuations, and Vs(k) is the 
power spectrum for the matter fluctuations, 5 = Sp/p, 
at the centre of the void. For any observer, each Vi will 
be normally distributed with mean zero and variance ob- 
tained by integration of the above spectrum: 

a ' 2 = /^ (fc) f 

= —p^jf Vs(k) ¥ , (5) 

where p c is the critical density, the FLRW scale factor 
a is set to unity today, and the integral is cut off at the 
scale of the Local Group, k ~ 1.0 Mpc -1 (the integral is 
only weakly sensitive to the precise cutoff value) . 

We obtain Ps(k) from the public Boltzmann code 
CAMB [26] by using parameters corresponding to an ef- 
fective open FLRW model that gives the same physics at 
recombination and evolution along the central worldline 
as the void model (see [2] for details). Importantly, in 
this approach we must assume that the primordial am- 
plitude at the centre of symmetry is the same as that 
at the LSS, which is directly constrained by the CMB. 
Although we do not have an explicit model for the for- 
mation of such a large void, we only consider growing 
mode void profiles in this work, so that at early times 
the void can be accurately described by linear theory on 
FLRW. Thus, our assumption states that we can con- 
sider the standard inflationary primordial spectrum to 
be superposed on the early linear void. Note also that 
Eq. (5) agrees with the dominant part of a calculation of 
the dipole using linearized general relativity [27]. 
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The observed dipole will then be the sum of induced 
and stochastic peculiar velocities. We can treat this as a 
vector sum, letting one component of the stochastic ve- 
locity, say Hi, lie along the direction of v T , with the other 
two components orthogonal. Then, the observed dipole 
will have magnitude v equal to the Euclidean norm of the 
vector (v T + V\, V2, V3). To derive the probability density 
function (pdf) for v, we recall that for a random variable 
X with pdf fx(x) (the lower-case argument denoting a 
particular value of the corresponding upper-case random 
variable), the pdfs for Y = X 2 and Z = \[X are given 
respectively by [28] 

Sv{y) = ±=. [fx(Vy) + fx(-Vv)} , (6) 

f z (z) = 2zf x (z 2 ). (7) 

With each normally distributed as above, applying (6) 
to each component of (v T + Vx, 1)2, U3) gives 

/(^)-(«) = -^e-^/^wf W 5 ) , 
\J2-nxo \ a J 

Adding the three squared components (which entails per- 
forming convolutions between their distributions) and 
applying (7) to the result yields the pdf for v, found to 
have the form of a "skewed Maxwellian:" 

f(v) = y?JL e -C-K)/(V) sinh(^) , v > 0. (8) 

V 7T (7V T V (7 / 

This can be integrated analytically to yield the cumula- 
tive density function for v in terms of error functions, a 
lengthy expression that we omit here. 

C. Choice of measured dipole 

Before using a measured dipole to constrain our posi- 
tion, careful thought should be given to which measured 
dipole to use. References [3, 14] use the dipole observed 
for the solar system's motion relative to the CMB sky 
(that is, the raw measured dipole after subtracting the 
motion of the detector relative to the Sun), v ~ 369 km/s. 
Meanwhile, Ref. [20] uses a w — 10~ 3 (v = 150 km/s) and 
Ref. [21] uses a w = 1.23 x 10~ 3 (v = 185 km/s), possibly 
due to confusion about the relationship between v and 
aio- Since we employ a linear treatment of fluctuations, 
we should ignore nonlinear contributions to the dipole, 
which are important only on very small scales. Also, we 
must consider that the entire Local Group, which ex- 
tends in space for roughly 1 Mpc, would be contained 
well within a void of the scale considered in this work. A 
second observer positioned in another galaxy within this 
group would therefore be expected to obtain the same 
limit as us as to how far from the void centre they could 
be. For these reasons, we choose the dipole measured 
for the Local Group itself, i>lg — 627 km/s [29], as the 
constraining value in our study. 
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FIG. 1: Maximum temperature anisotropy, \ AT\/T, vs. radial 
coordinate of the void centre for a selection of profiles from the 
chain (with no outer shells), with the grayscale level indicating 
the relative likelihood of the fit to CMB + SN data. Apart 
from where AT changes sign at r/r^ss — 0.8, the magnitude 
of the effect always exceeds 10 -3 . 



III. COMPUTATIONAL METHODS AND 
RESULTS 

A. Multiple voids 

As an accessory to the constraint in Sec. Ill C on hab- 
itable volume within a void, we first consider whether 
there could be other large voids between us and the LSS. 
If so, they could be expected to leave prominent scars 
on the CMB sky. This has previously been examined 
for voids that do not attempt to mimic the SN evidence 
for acceleration (see, e.g., [30-33]), and also for so-called 
"swiss cheese" models [34], but our analysis will take 
place within a slightly different framework. A spherical 
void would produce an axially symmetric anisotropy pat- 
tern on the sky, with a profile determined by the void's 
radial profile, and with an angular diameter dependent 
on the void's size and distance from us. The anisotropy 
would be strongest in the centre of the pattern, corre- 
sponding to photons that travel from the LSS through the 
void centre. We calculated the temperature anisotropy 
for photons passing radially through a void in Sec. II A, 
and we can reuse that technique here, but now placing 
the observer at distances from the void as great as the 
radius of last scattering, tlss 1 which is the maximum dis- 
tance at which a void could be visible. In the parlance 
of Sec. II A, the "outgoing" ray passes through the void, 
while we use the "incoming" ray to define the reference 
temperature T as a matter of convenience. 

While the voids examined in [2] required an outer over- 
dense shell to match CMB data, this would create a much 
larger effect than for just a central void. Therefore, for a 
conservative estimate of the effect of other voids, we dis- 
card the shell for this calculation, leaving only the cen- 
tral underdensity in an Einstein-de Sitter (EdS) back- 
ground. A plot of the resulting amplitude of the tern- 
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perature anisotropy, \AT\/T [calculated using Eqs. (1) 
and (2)] for a selection of void profiles from our Markov 
chain (described in Sec. Ill B) is shown in Fig. 1, with the 
grayscale level indicating the likelihood of the fit to CMB 
+ SN data. Up until r/r^ss — 0.2, the local dipole effect 
described in Sec. II A dominates, while past r/r^ss — 0.8, 
the void intersects the observer's LSS, resulting in the 
regular Sachs- Wolfe effect. Between these regions, a non- 
linear integrated Sachs- Wolfe (ISW), or Rees-Sciama, ef- 
fect contributes a \AT\/T with magnitude greater than 
10~ 3 . In comparison, the WMAP Cold Spot has roughly 
\AT\/T ~ 4 x 1CT 5 (see, e.g., [35]), so the presence of 
other large voids (i.e., of the kind required locally to fit 
the SN data) within our observable volume is ruled out 
by observations. In a recent study, no evidence was found 
for a void of radius ~ 200h~ 1 Mpc as a possible source 
for the Cold Spot itself [36]. 



B. Radial constraint 

We use COSMOMC [37] to generate Markov- Chain 
Monte-Carlo chains of void models as described in [2], 
distributed according to fits to the full CMB tempera- 
ture power spectrum (using data from WMAP7 [38, 39], 
ACBAR [40], BOOMERANG [41], CBI [42], and QUAD 
[43]) and SNe la (adopting the Union2 compilation [44], 
which covers the range z = 0.015-1.4). Each model is 
parametrized by its radial density profile 5{t,r) using a 
sum of polynomials representing a central underdensity 
and an overdense shell located at z > 1.5 (see [2] for ex- 
plicit expressions), and approaches an EdS background 
at large distances. These models consist of pure growing 
mode profiles, and hence are specified by a single radial 
profile. 

At a given proper distance from the centre (measured 
along a comoving slice), we evaluate the local dipole in- 
duced by each void model using Eq. (1), and then, using 
Eq. (8) , calculate the probability that the observed dipole 
v (which includes the stochastic peculiar velocity) will be 
less than ulg- The resulting value is checked by building 
a normalized histogram from samples of the distribution 
(8) and integrating up to v^q. Averaging the probabil- 
ities over the entire chain of models, and repeating this 
process for different physical distances, allows us to ob- 
tain confidence limits on how close to the centre of the 
void an observer must be in order to see a dipole less 
than vlg- 

We find that an observer must be within 80 Mpc of 
the centre of a void in order to observe a dipole less than 
i>lg with greater than 95% confidence. The loosening of 
the constraint on radial distance in comparison with past 
studies is caused by the three additional considerations 
described above: our broad exploration of the space of 
profiles; our choice of the Local Group dipole instead of 
a value based on smaller-scale motion to use as the em- 
pirical value; and the inclusion of the extra stochastic 
peculiar velocity component. If, for example, we use the 
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FIG. 2: Distribution of dipole magnitude v at radii 24 Mpc 
(green curves) and 80 Mpc (red curves) over the chain of 
void models, with random peculiar velocity (solid colours) 
and without (faint colours, apparently noisy due to the finite 
Markov chain sampling). At each distance, curves are scaled 
to the same height. The addition of peculiar velocity is seen 
to significantly widen both distributions. 

dipole from [20], namely 150 km/s, then the third con- 
sideration alone results in a radial constraint of 24 Mpc, 
pictured in green in Fig. 2, while the constraint using «lg 
is seen in red. Despite being significantly looser than pre- 
vious estimates, this new limit on radial distance using 
the CMB dipole anisotropy is still much tighter than that 
imposed by modern supernovae surveys [14, 15]. 



C. Volume constraint 

In addition to endeavouring to limit an observer's dis- 
tance from the centre, it is also worthwhile to ask what 
fraction of the total volume of a void one could fall within 
and still see a dipole less than i>lGj since this fraction 
takes into account the physical extent of the void itself. 
However, we must also consider the prior probability of 
finding ourselves within a void at all, and this depends on 
whether other voids of a suitable size can exist between 
our current position and the LSS. 

As detailed in Sec. Ill A, there can be a maximum of 
one void within our Hubble radius (namely, the central 
void which can fit the SN data). Therefore, we should 
actually seek the fraction of our entire Hubble volume 
we could reside in to give an appropriate dipole. If 
we imagine a process of placing an observer at a uni- 
formly random position, this fraction can be interpreted 
as the probability that the observer would end up seeing a 
dipole that agrees with i>lg- It is more sensible to picture 
this process occurring at an early time, so we evaluate the 
corresponding volumes on a spacelike slice corresponding 
to the time of last scattering. The process for obtaining 
confidence limits on the volume proceeds similarly to the 
radial calculation, described in Sec. IIIB. 

The result is that an observer must inhabit a fraction 
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FIG. 3: Fractional difference between luminosity distance in 
ACDM and in a void whose redshift-luminosity distance re- 
lation matches well to ACDM at the current epoch, plotted 
at different observation times. Both relations are recalculated 
at each observation time. The discrepancy exceeds 2% at ±1 
Gyr, and increases at later or earlier times. The worse fit at 
low redshift stems from our constraint that the void profile 
be smooth at the centre. 
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FIG. 4: We simulate the (x—z relationship for ACDM at var- 
ious times given by 2000 SNc uniformly distributed from 
z = 0.01-1.7, with a Gaussian dispersion characterized by 
<7 M = 0.15. We then calculate \ 2 values at those times for 
an effective "fit" of this relationship to that of a void profile, 
averaged over 30 fiducial data sets. There is a clear preference 
for the current epoch, with times more than ~1.5 Gyr into 
the past or future exhibiting severe deviations from ACDM. 



of the Hubble volume of order 10 -8 , again in order to 
have v < «lg a t 95% confidence. This is a clear violation 
of the Copernican principle, since observations force our 
position within the observable universe to have been fine- 
tuned to about one part in 100 million at early times. 
If we live in a void, our position becomes very special 
indeed. 



IV. TEMPORAL POSITION 

In light of the constraint on where we must live within 
a void, it is natural to ask if there might also be some re- 
strictions on when we live. It has often been pointed out 
that the SMC suffers from such temporal tuning, in that 
the densities of matter and dark energy are comparable 
today but begin to differ by many orders of magnitude 
in the relatively near past or future (see, e.g., [45]). Any 
void model has an inherent time dependence, since the 
radial density profile "flattens out" into the past while 
the central underdensity deepens and the outer shell be- 
comes even denser into the future. Thus, it can be ex- 
pected that observational relations should be drastically 
altered at different times. 

It is possible to investigate this by calculating the 
redshift-luminosity distance relation for a void as seen by 
an observer at different times, and comparing it to that 
of the SMC (otherwise known as ACDM) at those times. 
(We emphasize that relations in both void and ACDM 
models are time-dependent, so both must be recalculated 
for each observation time.) In Fig. 3, we examine a void 
profile that matches ACDM well at the current epoch, 
and find that the fractional difference between luminos- 
ity distances of the void and ACDM increases for times 



other than the present, as expected. 

Attempting to generate fiducial data can also lend in- 
sight into how observations in a void would change with 
time. To this end, we simulate 2000 SNe, uniformly dis- 
tributed in redshift from z — 0.01-1.7, by calculating the 
distance modulus /i for each one according to ACDM at a 
certain time and adding on a realistic Gaussian dispersion 
with cr^ = 0.15. We then find the x 2 of an effective "fit," 
treating the fi-z curve at that time from the void profile 
used in Fig. 3 as the null hypothesis to be compared with 
the simulated ACDM data. This process is repeated at 
different observation times, and the results are shown in 
Fig. 4. In addition to the best fit occurring at the cur- 
rent time, we see that observing the void more than ~1.5 
Gyr into the past or future results in A% 2 3> v / 2000 ~ 45, 
giving a redshift-luminosity distance relation that notice- 
ably deviates from the behaviour of ACDM at that time. 

Thus, a void that fits the relation of standard ACDM 
does so only temporarily. It would therefore be rather co- 
incidental to find ourselves living at a time when a void 
matches ACDM so closely. Of course, better or worse 
simulated data would yield a more or less obvious tem- 
poral fine-tuning. The important point, however, is that 
void models are not free of the temporal tuning that ap- 
pears to plague ACDM. 



V. CONCLUSIONS 

Our primary goal in this paper has been to obtain 
rigorous limits on the region an observer could inhabit 
within a cosmic void while detecting a CMB dipole less 
than the measured value. We have examined a wide 
range of void profiles for this purpose and clarified the 
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best choice for the measured dipole. It is conceivable that 
a peculiar velocity in the proper direction could cancel a 
large dipole induced by a void, and we have accounted 
for that possibility in a probabilistic fashion by including 
an extra velocity distribution for the observer, derived 
from the expected primordial perturbation spectrum. 

We found that an observer must be located within 
80 Mpc of the void centre, or within a fraction of the 
Hubble volume of order 10~ 8 , to see a dipole less than 
«lg at 95% confidence. This constraint is weaker than 
previous estimates, both as a radius and as a fraction 
of the void radius (voids in our chain have average radii 
of roughly 3 Gpc, giving an average constraint of ap- 
proximately 2.5% of the void radius, as opposed to the 
1% found in [14]). We also showed that there can be a 
maximum of one void large and deep enough to match 
the SN data within our Hubble volume by examining the 
anisotropy any other voids would create in the CMB. We 
additionally highlighted that the redshift-luminosity dis- 
tance relation of a void can resemble that of ACDM for 
only a brief epoch, and thus the time at which we observe 
a void is subject to some degree of fine-tuning. 

These results indicate tension with the Copernican 
principle, severe in space but also present in time, as an 
unavoidable consequence of void models to explain cos- 
mic acceleration. When coupled with the low probability 
of such a void forming by chance through the process that 
created the other matter perturbations [7] , it appears ex- 
ceedingly unlikely that the Universe would arrange itself 
in such a configuration. 

A natural question is whether departures from spheri- 
cal symmetry, i.e. moving beyond the LTB model, could 
alleviate the spatial tuning we have examined here. (This 
important question was addressed very recently in [46].) 
Because a void required to fit the SN data must be Gpc- 
scale, void models with departures from a spherical pro- 
file on the order of Gpc might substantially reduce the 



spatial tuning problem, simply because the central posi- 
tion would not be well defined in such models. However, 
nonlinear structures on Gpc scales (i.e. the features pro- 
ducing the departures from spherical symmetry) should 
produce an obvious nonlinear ISW imprint on the CMB, 
as we highlighted in Sec. Ill A. Therefore, while such 
models which depart significantly from spherical might 
lessen the tuning problem from the dipole, they would be 
expected to generate considerable power in higher multi- 
poles of the CMB. It appears inevitable that any depar- 
tures from spherical must be on scales much smaller than 
Gpc, in which case the spatial tuning problem persists. 

While these considerations may be seen as more philo- 
sophical than concrete, they complement recent thorough 
comparisons of void models with a wide variety of ob- 
servational data that find strong discrepancies between 
the two [2, 3] (however, see also [4]). As an example, 
in [2] alone it is found that while voids can be made to 
fit SN and CMB data, they possess Hubble rates, ages, 
radial BAO scales, Compton y-distortions, and local mat- 
ter fluctuation amplitudes that all come into conflict with 
observed values. There are many prospects as well for fu- 
ture constraints, including improved BAO measurements 
from surveys such as PAU [47] and BOSS [48], more pre- 
cise CMB spectra from Planck [49], and perhaps even 
other limits on spatial position from cosmic parallax and 
rcdshift drift [50, 51] . It is clear that void models continue 
to face serious obstacles to becoming successful alterna- 
tives to the SMC. 
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